Search for universal roughness distributions in a critical interface model 
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We study the probability distributions of interface roughness, sampled among successive equi- 
librium configurations of a single-interface model used for the description of Barkhausen noise in 
disordered magnets, in space dimensionalities d — 2 and 3. The influence of a self-regulating (de- 
magnetization) mechanism is investigated, and evidence is given to show that it is irrelevant, which 
implies that the model belongs to the Edwards- Wilkinson universality class. We attempt to fit our 
data to the class of roughness distributions associated to 1//" noise. Periodic, free, "window", and 
mixed boundary conditions are examined, with rather distinct results as regards quality of fits to 
1//" distributions. 

PACS numbers: 05.6.5. +b, 05.40.-a, 75.60.Ej, 05.70.Ln 



I. INTRODUCTION 

This paper deals with fluctuation properties of driven 
interfaces in random media. The subject has been the 
focus of much current interest (for reviews see, e.g., 
Refs. Q, Q). Special attention has been given to fea- 
tures at and close to the depinning transition, where a 
threshold is reached for the external driving force, above 
which the interface starts moving at a finite speed. In 
analogy with the well-established scaling theory of equi- 
librium critical phenomena, one usually searches for the 
underlying universality classes and their respective crit- 
ical indices, wherever such concepts are applicable. One 
example is the roughness exponent C which character- 
izes the disorder-averaged mean-square deviations of the 
interface about its mean height, at depinning []|. 

It has been shown very recently that the probabil- 
ity distribution functions (PDFs) of critical fiuctuations 
in seemingly disparate (both equilibrium and out-of- 
equilibrium]_sjst_ems display a remarkable degree of uni- 
versality [3, 4, 5, 6]. In the context of depinning phenom- 
ena, this indicates that one may gain additional insight 
into the physical mechanisms involved, by investigating 
the full roughness PDFs instead of concentrating on their 
lowest-order moments. Here we investigate the PDFs of 
interface roughness for a specific single-interface model 
which has been used in the description of Barkhausen 
noise [7, 8, 9, 10|, and is related to the quenched Edwards- 
Wilkinson universality class [ill Il2 . U^ ll4J . A pre- 
liminary investigation of this problem was reported in 
Ref.lifl 

Barkhausen "noise" (BN) is an intermittent phe- 
nomenon which reflects the dynamics of domain-wall mo- 
tion in the central part of the hysteresis cycle in ferromag- 
netic materials (see Ref . [ij for an up-to-date review) . A 
sample placed in a time-varying external magnetic field 
undergoes sudden microscopic realignments of groups of 
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magnetic moments, parallel to the field. For suitably slow 
driving rates, such domain- wall motions, or "avalanches", 
are well separated and can be easily individualized. The 
accompanying changes of magnetic fiux are usually de- 
tected by wrapping a coil around the sample and measur- 
ing the voltage pulses thus induced across the coil. The 
integral of the voltage amplitude of a given pulse over 
time is proportional to the change in sample magnetiza- 
tion, thus giving a measure of the number of spins over- 
turned in that particular event, or "avalanche size". Mod- 
ern experimental techniques allow direct observation, in 
ultra-thin films, of the domain- wall motion characteristic 
of BN, via the magneto-optical Kerr effect (ln.ll7l|. 



It has been proposed that BN is an illustration of "self- 
organized criticality" P, Ha, llil l2(l |. in the sense that 
a broad distribution of scales (i.e. avalanche sizes) is 
found within a wide range of variation of the external 
parameter, namely the applied magnetic field, without 
any fine-tuning. Accordingly, the interface model stud- 
ied here incorporates a self-regulating mechanism in the 
form of a demagnetizing term (see below). In the con- 
text of interface depinning models, the question arises 
of whether this is a relevant perturbation, i.e., whether 
self-organized depinning phenomena belong to the same 
universality class as their counterparts which do not in- 
corporate such mechanisms. 



In what follows, we first recall pertinent aspects of the 
interface model used here, and of our calculational meth- 
ods, as well as the connections between roughness dis- 
tributions and 1//" noise. Next, we exhibit numerical 
data for roughness distributions, generated by our simu- 
lations. We examine the infiuence of the self-regulating 
mechanism, and investigate the effect of assorted bound- 
ary conditions, both on our results and on the class of 
1//" noise distributions to which they are compared. Fi- 
nally, we discuss our findings with regard to the relevant 
universality classes. 



II. MODEL AND CALCULATIONAL METHOD 

The single-interface model used here was introduced in 
Ref. Q for the description of BN. We consider the adia- 
batic limit of a very slow driving rate, thus avalanches 
are considered to be instantaneous (occurring at a fixed 
value of the external field) . 

Simulations are performed on an Lx x Ly x co ge- 
ometry, with the interface motion set along the infinite 
direction. The interface at time t is described by its 
height hi = h{x,y,t), where {x,y) is the projection of 
site i over the cross-section. No overhangs are allowed, 
so h{x,y,t) is single- valued. We consider mainly Ly = 1 
(system dimensionality d = 2, interface dimensionality 



d' 



1 = 1), and Lx — Ly {d — 3, d' 



For 



reasons to be explained below, we will use the following 
sets of boundary conditions: periodic (PBC), so every 
site has two neighbors for d = 2 and four for d — 2t\ free 
(FBC), meaning that the interface is horizontal at the 
edges (dh/dn — 0, where n = 2? or y is the normal in the 
cross-section plane), and mixed (MBC), i. e., periodic 
alongx and free along y. These latter were employed in 
Ref. |l3, to reproduce the physical picture of films with 
varying thickness. We also considered an alternative im- 
plementation of FBC, namely window boundary condi- 
tions (WBC), to be described in Section ITV CI 

Each element i of the interface experiences a force given 
by: 



h '^u{x,y,hi) + k^ [h. 



j 



,W 



H, 



where 



He = H -f]M . 



(1) 



(2) 



The first term on the RHS of Eq. 1^ is chosen randomly, 
for each lattice site fl = {x,y,hi), from a Gaussian dis- 
tribution of zero mean and standard deviation R, and 
represents quenched disorder. Large negative values of u 
lead to local interface pinning. The second term (where 
the force constant k is taken as the unit for /) corre- 
sponds to elastic nearest-neighbor coupling (surface ten- 
sion); £j{i) is the position of the j-th nearest neighbor of 
site i. For MBC, sites at y = 1 and y = Ly have only 
three neighbors on the xy plane (except in the monolayer 
case Ly = 1 which is the two-dimensional limit, where all 
interface sites have two neighbors). The last term is the 
effective driving force, resulting from the applied uniform 
external field H and a demagnetizing field which is taken 
to be proportional to M = (l/LxLy) J2i=i " ^i^ the mag- 
netization (per site) of the previously fiipped spins for a 
lattice of transverse area LxLy. For actual magnetic sam- 
ples, the demagnetizing field is not necessarily uniform 
along the sample; even when it is (e.g. for a uniformly 
magnetized ellipsoid) , rj would depend on the system's 
aspect ratio J2y. Therefore, our approach amounts to 
a simplification, which is nevertheless expected to cap- 
ture the essential aspects of the problem [2|. Here we 



use R = 5.0, k = 1, rj = 0.05, values for which fairly 
broad distributions of avalanche sizes and roughness are 
obtained [3, 1^, Il3| ■ We also consider the effects of taking 
77 = 0, i.e., the non-self-organizing limit. 

We start the simulation with a fiat wall. All spins 
above it are unfiipped. The force fi is calculated for 
each unfiipped site along the interface, and each spin at 
a site with fi > fiips, causing the interface to move 
up one step. The magnetization is updated, and this 
process continues, with as many sweeps of the whole lat- 
tice as necessary, until fi < for all sites, when the 
interface comes to a halt. The external field is then in- 
creased by the minimum amount needed to bring the 
most weakly pinned element to motion. The avalanche 
size corresponds to the number of spins fiipped between 
two consecutive interface stops. 

On account of the demagnetization term, the effective 
field He at first rises Hnearly with the applied field H, 
and then, upon further increase in H, saturates (apart 
from small fiuctuations) at a value rather close to the 
critical external field for the corresponding model without 
demagnetization 0,01 ■ The saturation Hf. depends on i?, 
k and 77 {not noticeably on Lx, Ly) ja, H^, and can be 
found from small-lattice simulations. It takes 10^ — 10'^ 
avalanches for a steady-state regime to be reached, as 
measured by the stabilization of He against H. 



III. ROUGHNESS DISTRIBUTIONS AND 1//" 
NOISE 

We have generated histograms of occurrence of inter- 
face roughness, to be examined in the context of universal 
fiuctuation distributions |^ H, IM lOl - We have used only 
steady-state data, i.e., after the stabilization of He of 
Eq. (El against external field H. This is the regime in 
which the system is self-regulated at the edge of critical- 
ity 0,1^' As the model is supposed to mimic the data 
acquisition regime for BN during which the external field 
grows linearly in time 0) la S Efl Hal ) the value of H is 
a measure of "time". 

At the end of each avalanche, we measured the rough- 
ness -0)2 of the instantaneous interface configuration at 
time t, as the (position-averaged) square width of the 
interface height |ftl22|: 



hj^Ljj 



W2{t) ^ [LxLy)-^ J2 {h{t)~h{t)Y 



(3) 



where h{t) is the average interface height at t. As the 
avalanches progress, one gets a sampling of successive 
equilibrium configurations; the ensemble of such config- 
urations yields a distribution of the relative frequency of 
occurrence of W2 ■ Here we usually considered ensembles 
of 5 X 10^ events (on e and a half orders of magnitude 
larger than in Ref. IkI) , so we ended up with rather clean 
distributions. This was essential, in order to resolve am- 
biguities left over from our previous results [l(l |. 



The width distributions for correlated systems at crit- 
icality may be put into a scaHng form [fll Ift 122. 12^ . 



$(z) = (W2)^(W2) , 2 = ^2/(^2), 



(4) 



where angular brackets stand for averages over the en- 
semble of successive interface configurations, and the 
size dependence appears only through the average width 
(102). By running simulations with 0(10^) events, and 
400 < Lx < 1200 for rf = 2 {Ly = 1), m < L^ = Ly < 80 
for d = 3 ilO|, we ascertained that Eq. Q indeed holds, 
i.e., finite-size effects are not detectable in any significant 
way as far as the scaling functions '^(z) are concerned. 
The finite-size scaling of the first moment gives the rough- 
ness exponent [l|: 



{W2{L))^L^^ , 



(5) 



In the context of critical fiuctuation phenomena, it is 
known that boundary conditions have a non-trivial effect 
on scaling functions, as infinite-range critical correlations 
are sensitive to the boundaries of the system [a, lu |2^ 
124 12^ . This is the motivation for use of the assorted 
boundary conditions defined in SecHll 

We have compared our results against the family 
of roughness distributions for 1//" noise, described in 
Refs. 0, |2^ As explained there, such distributions are 
derived under the assumption that the Fourier modes 
into which the interface is decomposed are uncorrelated 
(generaHzed Gaussian approximation J23|), and with am- 
plitudes such that the frequency dependence of the power 
spectrum is purely 1//" [a|. This is the simplest starting 
point from which one may expect non-trivial results (the 
trivial ones corresponding to the case in which the real- 
space fluctuations are themselves uncorrelated, implying 
a =1/2). 



IV. RESULTS 

A. Influence of self-regulating term 

We first investigated what could be learned about the 
relevance of the self-regulating term, as regards rough- 
ness distributions. In order to do so, we determined the 
approximate critical value H^ of the internal field He of 
Eq. l(2l, by starting a simulation with rj ^ and wait- 
ing for He to stabilize. At that point, we set 77 = and 
repeatedly varied H in the interval (xiJg,i/g), a; < 1, 
according to the procedure dehneated in SeclH Though 
the interval of variation of H did affect the size distri- 
bution of avalanches, as this is what characterizes the 
proximity of the depinning point [tI, 1^, no change was 
apparent in the roughness data when comparing results, 
e.g., for X = 0.95 and x = 0.9. For the simulations de- 
scribed in the remainder of this subsection, we used the 
latter value. In all cases studied, namely, d = 2 PBC 
and d = 3 with both MBC and PBC, the infiuence of the 
demagnetization term on the roughness PDFs is rather 
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Figure 1: (a) Scaled probability distributions $(2:) in d = 3 
with MBC, for z defined in Eq. ^. Data for L = 40, 5 x 10^ 
configurations. Full line: demagnetization factor 77 = 0.05; 
dashed line: 77 = 0. (b) Scaling function difference against z. 



small, but systematic. This is illustrated in Fig. Q] for 
d = 3 with MBC, the case for which the deviations be- 
tween the 77 ^ and 77 = sets of data are the largest in 
magnitude. One sees that neglecting the demagnetizing 
term causes a small leftward shift of the scaHng curve. 



As we will see in Section lIV BL the changes it causes to 
the fits of our distributions to the analytical 1//" curves 
are of the order of systematic imprecisions characteristic 
of the fitting procedure. Nevertheless, it is instructive to 
seek the physical origins of such effect. This is done by 
direct inspection of the unsealed PDFs. In Fig. 12 it is 
apparent that, for 7; = the high-end tail of P{w2) is 
sHghtly fatter than for 77 7^ 0, at the expense of a small 
amount of depletion around the most probable value of 
u'2. Accordingly, the average {W2) is higher by ~ 8% 
in the former case than in the latter (the fractional dif- 
ference between averages is the same also for d = 2 and 
d = 3 PBC). Such a trend can be understood by recalling 
that the 77 = data have been collected just helow the 
depinning transition, i.e., still within the regime where 
pinning forces are dominant. Thus the interface mostly 
meanders about, in order to comply with local energy 
minimization requests. The complement of this picture 
is that, for H > He the interface moves with finite speed, 
more or less ignoring local randomness configurations, 
and becoming smoother the farther one is above the crit- 
ical point. In short, for a given lattice size the average 
interface roughness decreases monotonically as the ex- 
ternal field (driving force) is increased across its critical 
value. 
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Figure 2: (a) Probability distributions P{w2) in d = 3 with 
MBC. Data for L = 40, 5 x lO'^ configurations. Full line: 
demagnetization factor "q = 0.05; dashed line: rj = Q. (b) 
Probability distribution difference against W2. Extent of hor- 
izontal axis corresponds to the same interval of z-variation in 

Fig.m 



The interpretation of the small differences between 77 = 
and 77 ^ distributions is then as follows: (i) because of 
the way in which data for the former were collected here, 
they represent a system just below He, for which interface 
roughness is sHghtly larger than at the critical point; and 
(ii) the closeness of 77 = data to those for ry ^ 0, and 
the way in which both sets of data differ, strongly suggest 
that behavior at the critical point of the 77 = system is 
the same as that of the 77 7^ (self-regulated) case. We 
conclude that the self-regulating term is irrelevant, as far 
as critical roughness distributions are concerned. 



B. PBC, d=2and3 

Analytical expressions for the 1//" distributions with 
PBC are either given in Ref.@ (d = 2), or can be derived 
straightforwardly from Refs. |a, |2^ {d = 3). In the lat- 
ter case, the use of exact identities for two-dimensional 
lattice sums pfl| speeds up calculations considerably. 
Estimates of the exponent C, of Eq. Q, from power- 
law fits of simulational data with 0(10^) events, and 
400 <L.^< 1200 for d = 2, 30 < L < 80 for d = 3_^ive 
C(d = 2, PBC) = 1.24(1), C(rf = 3, PBC) = 0.71(1) |l3|. 

Consideration of the scaHng properties of height-height 
correlation functions and their Fourier transforms then 
suggests [221, for the generalized Gaussian case of inde- 



Figure 3: x^ per degree of freedom (x^d.o.f.) for fits of simula- 
tion data with PBC to analytical forms of 1//" distributions, 
against a. Triangles: d = 2, Lx = 400; squares, d = 3, 
L = 40. 



pendent Fourier modes, that 

a^d' + 2C {d' 



(6) 



which would imply a = 3.48(2) {d = 2), 3.42(2) {d = 3). 

Such predictions can be quantitatively checked by es- 
timating the values of x^ per degree of freedom (x^^ ^ f ) 
from fits of our simulation results to the analytical dis- 
tributions. Since, even with 5 x 10^ samples, the simu- 
lational data eventually get frayed at the top end, given 
the long forward tails characteristic of all systems stud- 
ied here, our fits used only data for which ^(z) > 10~^. 
This turned out not to be a drastic restriction, as we 
were left typically with at least 100 — 200 points to fit 
in each case. Assuming the uncertainty in the value of 
a that best fits our data to be given by requiring that 
X^d o f ^tay within 150% of its minimum, we quote from 
the' data shown in Fig.El a = 3.60(13) (d = 2); 3.52(6) 
(d = 3). The agreement with the above predictions is 
satisfactory, though slight discrepancies remain. A vi- 
sual check of the goodness-of-fit for each case is given in 
Figs. Hand El 

Fitting 77 = data to the closed-form distributions 
produces curves whose minima of x^d o f ^^^ essentially 
the same as in Fig. El and slightly shifted rightwards. 
Using the same criteria as above for the estimation of 
error bars, we have, for 77 = 0: a = 3.64(16) {d = 2); 
3.59(5) (d = 3). 

Detailed discussion, and pertinent comparisons with 
data from Ref. [22, will be deferred to Section IVl 
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Figure 4: Scaled probability distribution $(z) in rf = 2 (PBC), 
for z defined in Eq. Q, from 5 x 10^ configurations. Squares: 
simulation data {L = 400). Full line is roughness distribution 
for 1//" noise given in Ref. [y, with a — 3.60. 



C. FBC and WBC, d = 2 and 3 

We have generated roughness data in both d — 2 and 3 
with FBC. Our initial implementation of FBC, used also 
in Ref . llOl aims at a literal reproduction of the constraint 
that the interface must be horizontal at the edges. Thus, 
e.g. for d = 2, "ghost" sites are added at x = 0, x = Lx + 
1, whose heights are always adjusted to be respectively 
h{0,t) = h{l,t), h{Lx + l,i) = h{La,,t). This way, the 
edge sites at a: = 1 and L^ experience no elastic pull (see 
the second term on the right-hand side of Eq. (QJ) from 
their ghost neighbors outside the sample. 

Similarly to the PBC cases, estimates of the exponent 
C of Eq. were extracted from power-law fits of simu- 
lational data with 0(10^) events, and 400 < L^ < 1000 
for d = 2, 30 < L < 80 for d = 3. The results are 
C{d = 2, FBC) = 1.28(2), ({d = 3, FBC) = 0.89(1). 
While the former value might be construed as not incon- 
sistent with PBC and FBC giving the same universality 
class for d = 2, the same picture cannot hold for d = 3. 
Though it is known 5, 6, 22, 24,^3 that boundary con- 
ditions do have significant influence on scaling functions 
of critical systems, they are not generally expected to 
change the values of critical exponents. 

In order to discuss the roughness PDFs, we first recall 
the effect of FBC on 1//" distributions. The generating 
function G(s) = J dw2 P{w2) e~*'"^ has the general form 
for PBC Ulil 



Gpis) 







Figure 5: Scaled probability distribution $(z) in d = 3 with 
PBC, for z defined in Eq. JIJ, from 5 x 10 configurations. 
Squares: simulation data {L = 40). Full line is roughness 
distribution for 1/f" noise, with a — 3.52. 
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(7) 



where n is a lattice vector in d — 1 dimensions with in- 
teger coordinates. Because all n are counted, the square 
root disappears due to the (at least) twofold degeneracy. 
Requiring that the interface be horizontal at the edges 
implies that the Fourier representation of h{t) includes 
only cosines. The corresponding Gf{s) has the degener- 
acy of its singularities cut in half, compared to PBC. 

In d = 2, this means that the single poles found for 
PBC turn into square-root singularities. Evaluation of 
P{w2), as the inverse Laplace transform of Gf{s), thus 
necessitates a direct approach, since the residue theorem 
is inapplicable. This has been accomplished in Ref. |23, 
from which the relevant expressions were extracted in or- 
der to attempt a minimization of x^^ ^ f against a, similar 
to that of Section ITVBI With ({d = 2, FBC) as above, 
one would expect a good fit for a ~ 3.5 — 3.6. Instead, 
X^d o f ^^^ ^ minimum value ~ 4 x 10~^ at a = 2.96, and 
increases monotonically to reach ~ 4 x 10~^ at a = 3.5. 
This is clearly at variance with correspondings results for 
the PBC case. 

We then decided to generate data using window bound- 
ary conditions (WBC) ja, 123; which are generally ac- 
cepted as an alternative way to simulate free edges. Ac- 
cordingly, in d = 2 we imposed global PBC on a system 
of overall length L^, and measured the local roughness 
within each of Uw adjacent windows of length Lx/n^. 



With n^ ^ 1, it is plausible to assume that the resulting 
PDFs are independent of the boundary conditions estab- 
lished at X — 0, Lx. In order to guarantee statistical 
independence, one should in principle use widely sepa- 
rated windows. However, the use of nonoverlapping, but 
neighboring, windows instead appears to introduce no 
measurable errors on the resulting PDFs |a|. We fixed 
riw = 10, and initially measured C via Eq. ^, from a 
sequence of simulations with 0(10^) events (i.e. individ- 
ual avalanches, thus the total number of roughness sam- 
ples is larger by a factor of Uu,), and 400 < L^ < 1200, 
which gave ({d = 2,WBC) = 1.21(2). Though this 
differs by 3.5 standard deviations from the value com- 
ing from FBC, it is just consistent, at the margin, with 
C(d = 2,PBC) = 1.24(1) found above. 

Direct examination of scaled PDFs results in the fol- 
lowing observations. First, in Figure El one can see that 
the PDFs in d = 2 for FBC and WBC are unmistakably 
distinct. Furthermore, fits of FBC data to the analytical 
expressions derived in Ref . [23 have been found to be gen- 
erally of low quality. As mentioned above, the best fit of 
FBC data is for the a = 2.96 curve, shown in the Figure 
as a dashed line, and corresponds to x^d of — 4 x 10~^. 
Though this average deviation is of the same order as 
that for the best case with PBC (recall Fig. EJ, compar- 
ison to Fig. 31 shows that, while for PBC discrepancies 
are concentrated close to the narrow peak (thus they can 
be at least partially ascribed to binning effects) , here one 
has a rather widespread disagreement in shape. 

On the other hand, WBC data can be much more 
closely fitted by the analytical expressions, as shown both 
in the inset of Fig. JHl where x^d o f exhibits a minimum 
value ~ 7 X 10"'* at a = 3.85, and directly in the main 
Figure, by the superposition of the a — 3.85 curve onto 
the corresponding numerical data. 

In summary, an analytical form derived from assum- 
ing an interface whose Fourier representation has only 
cosines (i.e. is horizontal at the edges) has provided 
a very good fit to numerical data generated by impos- 
ing WBC. Though this appears contradictory, the same 
procedure has been successfully accompHshed in Ref. [23, 
with regard to both experimental and simulational data. 

Still, an important question remains, since the op- 
timum a = 3.85(5) (error bars estimated as in Sec- 
tion liVB)l impHes C = 1.43(3) via Eq. ©. This is signifi- 
cantly distinct from all three estimates thus far obtained 
for d = 2, which average to 1.25(5). We shall defer the 
discussion of this point to Section Fvi 

Turning now to d = 3, all poles of Gp{s) have even 
degeneracy. A straightforward adaptation for FBC is as 
follows. Recalling that the lattice sums J^ |n|~" which 
crop up in the calculation of (^2) [E M l2^ must be 
halved, this impHes a rescahng of the variable s, so for- 
mally one can write J^: 




Gf{s) = jGp{2.s 



(8) 



Fitting our d = 3 FBC data to analytical distribution 
functions, obtained with help of Eq. IJHJ, turns out to 



Figure 6: Scaled probability distribution $(2:) in d = 2, for 
z defined in Eq. (UJ. Points are simulation data. Crosses: 
L = 400, FBC, 5 X lO'^ configurations. Squares: L = 400, 
WBC, lO'^ avalanches, n„ — 10 windows. Full line is rough- 
ness distribution for 1/f" noise (see Ref. |25t), with a — 3.85. 
Dashed line : roughness distribution for a = 2.96 (see text). 
Inset: x d.o.f. against a, for fits of WBC simulation data 
against l/f" distributions, showing a minimum at q = 3.85. 



give similar results to the d = 2 case. The above-quoted 
value C = 0.89(1), from the finite-size scaling of (^2), 
together with Eq. Q, would suggest a = 3.78(2). How- 
ever, x^d o f against a has a single minimum (~ 10"'^) at 
a = 3.18(8) (error bars estimated as in Section lTV B|l and 
increases monotonicaly, reaching ~2xl0~^atQ; = 3.78. 
We again resorted to WBC. Imposing PBC at the edges 
of a system with L x L cross-section, we measured lo- 
cal roughness within each of n^ non-overlapping, ad- 
jacent, square windows of linear dimension L/^/n^ (or 
the largest integer contained in it). We took n„ — 16, 
and initially measured C from a sequence of simulations 
with 0(10^) events, and 30 < L < 80, which gave 
(^{d = 3, WBC) = 0.75(2). The discrepancy between 
this and the value 0.89(1) coming from FBC is rather 
more severe than the corresponding case for d = 2. On 
the other hand, the present estimate is close to the val- 
ues of C(d = 3, PBC) found above, namely 0.71(1) from 
Eq.O and 0.76(3) from optimization of fits against 1//" 
distributions plus Eq. El 

Again, we investigated the roughness PDFs generated 
with WBC. Similarly to the d — 2 case, they differ 
markedly from the ones obtained with FBC, as shown 
in Fig. . This time, fits against the analytical expres- 
sions given through Eq. l|S|l exhibit a deep, well-defined 
minimum of x^d of sX a = 3.76(5) (see inset in the Fig- 







Figure 7: Scaled probability distribution $(2:) in d = 3, for 
z defined in Eq. ^^. Points are simulation data. Crosses: 
L = 40, FBC, 5 X 10^ configurations. Squares: L = 40, WBC, 
3 X lO'^ avalanches, n^, = 16 windows. Lines are roughness 
distributions for 1//™ noise (see Eq. lO), with a = 3.76 (full) 
and 3.18 (dashed). Inset: X^d.o.f. against q, for fits of WBC 
simulation data against 1//" distributions, showing a mini- 
mum at Of = 3.76. 



ure), in very good agreement with a — 3.78(2) predicted 
from finite-size scaling of (^2) data for FBC, together 
with Eq. (jHJ. However, for reasons to be explained at 
length in Section EI we believe this coincidence to be 
accidental. 



D. MBC, d = 3 

We started by studying systems with a square cross- 
section, imposing PBC along x and FBC, as defined at 
the beginning of Section ITV CA along y. 

Estimates of the exponent ( of Eq. ^ were again 
extracted from power-law fits of simulational data with 
0(10^) events, and 30 < L < 80 for d = 3 MBC, with 
the result C{d = 3, MBC) = 0.87(1) [l3|. 

The Fourier representation of h{t) with MBC can be 
put in the form: 



H^,y) = 51' 



^2iri(m2;+f j,)/L 



(9) 



where m^n — 0, ±1, ±2, . . . , {m, n) ^ (0, 0), and C-m.n — 
(^*n,n> Cm-n = c^.n- Thus & global rcscaling such as that 
of Eq. (IHJ is not possible. On the other hand, starting 
from Eq. IJHJ, an analysis similar to that of R.efs. Im. l2^ 



suggests a generating function: 

Gm{s) = n fi 



(4m2 + n2)"/2 



-1/2 



(10) 



again with (m, n) ^ (0, 0). The double sum ^^ „(4m-^ + 
^2-j-a/2^ which appears in the subsequent expression for 
(W2), corresponds to Q(l, 0, 4) of Ref.SQ and can be eas- 
ily evaluated. 

We performed fits of simulational data to the closed- 
form PDFs calculated as above. While Eq. IJHJ, with C, — 
0.87(1), gives a = 3.74(2), x^d o f ^^^ ^ minimum ~ 2 x 
lO^'^ at a = 3.36(10). The overall quality of fits is slightly 
worse than for d = 3 FBC (refer to Fig.dl. 

In order to investigate WBC, we took rectangular sys- 
tems with dimensions L^ and Ly = 4Lx with full PBC 
(we denote this setup as mixed window boundary con- 
ditions (MWBC)) and calculated local roughness distri- 
butions within Uw = 4: square windows of Lx x Lx sites 
each, side by side along the y axis. Scaling of the first 
moment of the distribution, Eq. @, with 30 < L^; < 80, 
gave C = 0.74(1). 

Again, the roughness PDF thus obtained was markedly 
distinct from that with MBC. In addition, fits to the an- 
alytical expressions derived from Eq. Ijl0|l were consid- 
erably worse than those of MBC data, with a minimum 

X^d.o.f - 1 ^ 10"^ at a = 4.1. 

The results are displayed in Fig.|Hl where it can be seen 
that even the best-fitting analytical PDF fails to provide 
a good match to the MWBC data (except for the initial, 
rather steep, ascent close to z = 0). 



V. DISCUSSION AND CONCLUSIONS 

We begin our discussion by recalling from Ref. 'IQ and 
Sec. IIVBI that, for the model considered here with PBC, 
the finite-size scaling of the first moment of the distri- 
bution gives C{d = 2, PBC) = 1.24(1), ({d = 3, PBC) = 
0.71(1). Both compare well with the usually accepted 
values for the quenched Edwards- Wilkinson (EW) uni- 
versality class m 113, HE 13, respectively C — 1-25 
{d = 2) and C, — 0.75 {d = 3). Furthermore, consid- 
eration of the full distributions points the same way: our 
simulational data displayed in Figs. 21 and El match very 
well those in Figure 2 of Ref. 2,2; which concern the EW 
model. The agreement with EW behavior is consistent 
with our results of Sec. IIV Al regarding the independence 
of scaled roughness distributions on the demagnetizing 
term. Indeed, the quenched EW equation can be written 
as ilJI 



i9/i(x,t) 
dt 



u(x,/i) + aV^/i(x) + / 



(11) 



where u represents quenched disorder and / is the exter- 
nal driving force. This has a one-to-one correspondence 
with Eq. (QJ , except that in that Equation we allowed for 
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Figure 8: Scaled probability distribution <l?(z) in d = 3, for z 
defined in Eq. (Ujl. Crosses: simulation data {L = 40, MBC, 
5 X 10 configurations). Squares: simulation data (Lx = 40, 
L:, = 160, MWBC (see text), 3 x 10** avalanches, n^^ = 4 
windows). Lines are roughness distributions for 1//™ noise 
(see Eq. ^), with a = 4.1 (full) and 3.36 (dashed). Inset: 
X d.o.f. against a, for fits of MBC simulation data against 
1//" distributions, showing a minimum at a = 3.36. 



the self-regulating, demagnetizing, term. Having shown 
that such mechanism is irrelevant as far as scaled rough- 
ness distributions are concerned, it becomes tenable to 
assume that, overall, our model belongs to the EW uni- 
versality class. 

Still for PBC, the connection between the exponents 
a and C,, predicted [III in Eq. JSJ, is verified within rea- 
sonable error bars. 

Turning to different sets of boundary conditions, we 
first point out that small differences in implementation 
of FBC (namely, "literal" FBC, i.e. horizontal interface at 
the edges, versus WBC) significantly alter the roughness 
PDFs. The question then arises of which, if any, of these 
implementations is the "right" one. 

We investigate this by referring to results derived 
through a "proven" method, i.e. finite-size scaling of the 
first moment of the distribution. Examination of the cor- 
responding column of Table Ql strongly suggests that, in 
both d = 2 and 3, WBC (including WMBC) preserves 
universality with PBC, while FBC does not (though in 
d = 2 FBC does not perform very badly). Accepting such 
preservation as a basic tenet, we conclude that FBC as 
implemented induces strong distortions in the scaling be- 
havior of interface roughness. In this context, the good 
agreement in d = 3 between the optimum a for fits of 
WBC data to the analytical forms, and that coming from 



Table I: Estimates of roughness exponent C, for different di- 
mensionalities and boundary conditions (BC). C^^^: finite- 
size scaling of first moment of distribution, Eq. (O. ("fit: 



from best-fitting 1//° 



value of X 



for C 



distribution and Eq. JSJ. 

- Cfit- 



X d.o.f. (min) 



d = 2 PBC 
d = 2 FBC 
d = 2 WBC 
d = 3 PBC 
d = 3 FBC 
d = 3 WBC 
d = 3 MBC 
d = 3 MWBC 



1.24(1 
1.28(2 
1.21(2 
0.71(1 
0.89(1 
0.75(2 
0.87(1 
0.74(1 



Cfit 

1.30(8) 
0.98(7) 
1.42(3) 
0.76(3) 
0.59(4) 
0.88(1) 
0.68(5) 
1.05(10) 



X d.o.f. (min) 

6 X 10"* 
4 X 10"^ 

7 X 10-* 
3 X 10"* 

1 X 10"^ 

8 X 10"^ 

2 X 10"^ 
1 X 10-^ 



finite-size scaling of FBC data via Eqs. Q and (EJ, must 
be regarded as fortuitous. 

Thus, we discard FBC, as well as MBC, for the re- 
maining of the present discussion. One must note, how- 
ever, that use of MBC (i.e. partial FBC) provides a sen- 
sible representation of the physical setup found in thin 
films, as well as reproducing well-known results (concern- 
ing scaling behavior of avalanche sizes) at both ends of 
the crossover between d = 2 and 3 [T^l . 

Returning to roughness scaling, we see in Tabled that 
the fair agreement between C^^^ and Cfit, found for PBC 
in d = 2 and 3, is absent in the remaining cases under con- 
sideration, i.e. d = 2 WBC, d = 3 WBC, d = 3 MWBC. 
One might ask whether finite-size effects (though widely 
believed to vanish already for small lattices P.lfl. l22 . l23 |) 
still have a nonnegligible quantitative effect on the scaled 
roughness PDFs found here, so as to distort our fits to 
the analytical distributions. We present data to show 
that this is not the case. 

In Figure El we compare L = 40 and i = 80 PDFs, 
for d = 3 WBC. Contrary to the systematic trend exhib- 
ited in Fig. n (for comparison between 77 = and 7^ 
distributions), here the difference A$(z) is rather small 
and essentially random, arising because of fiuctuations 
in statistics, coupled with binning effects. An apparently 
systematic effect shows up only for the narrow range close 
to z = where both PDFs have a steep slope. That, 
however, involves only of order 5 — 10 points, with a con- 
sequently reduced effect on the overall statistics. The 
corresponding curves y^^ ^ j against a are nearly indis- 
tinguishable; with L = 80 data, the minimum of x^^j ^ f is 
9 X 10"^ at a = 3.76(4), virtually identical to the L'= 40 
result shown in Fig [3 (see also Table |lj|. For d = 2 WBC 
and d = 3 MWBC, the overall picture is the same. There- 
fore, finite-size effects on the numerically-obtained PDFs 
are not a likely source for the disagreements found. 

We note also that, when considering 1//" distribu- 
tions, there is no apparent reason why Eq. lO should not 
hold for boundary conditions other than PBC, as that 
Equation was derived for generalized Gaussian distribu- 
tions I22I with the only assumption that the large-scale 
behavior is determined by a single observable. 




Figure 9: (a) Scaled probability distributions $(2:) in d = 3 
with WBC, for z defined in Eq. (gjl. Squares: L = 40. Full 
line: L = 80. In both cases, 10 avalanches, n„ — 16 windows, 
(b) Scaling function difference against z. 



a long-range elastic term, \ dx\ {h{x) — h{xi)) /(x—xi)'^, 
instead of the local term, V^/i(x), present here. While 
in that case an 1//" distribution gives good fits to the 
numerically-generated roughness PDF with WBC, this 
does not necessarily imply that a similar quality of fit 
can be found for the present EW problem with WBC. 
In this connection, one might ask how far the indepen- 
dent Fourier mode assumption, basic in the derivation of 
1//" PDFs, is affected by such details. One sees that the 
long-range term contributes qualitatively in the same di- 
rection as PBC, i.e. by imposing additional constraints 
on interface roughness (when compared, respectively, to 
short-range interactions and WBC). 

A plausible scenario then emerges, in which the ampli- 
tude of corrections to the representation of an interface 
roughness PDF by an 1//" distribution would depend on 
how much that interface is constrained, either by bound- 
ary conditions or by elastic terms in the equation of mo- 
tion. Lessening of such constraints would imply an in- 
crease in the correction amplitudes. However, at present 
we do not see a way to quantify and test these remarks. 

Clearly, more work is needed in order to clarify the 
connection between 1//" distributions and generalized 
depinning transitions. 



We are thus left with a single point to analyze, namely 
the overall adequacy of 1//" distributions to describe the 
problem at hand. The following comments are in order: 

(1) already for PBC, the study of generaHzed depin- 
ning problems shows that small but systematic discrep- 
ancies remain between numerical data and 1//" PDFs, 
whose origins can be traced to higher cumulants of the 
correlation functions J22|- Thus, in this sense the 1//" 
distributions are not expected to be a perfect fit, even 
for PBC. 

(2) In R ef ■ 1271 the equation of motion for h{x) contains 



Acknowledgments 

The author thanks Tibor Antal and Zoltan Racz for 
their advice on numerical evaluation of the closed-form 
PDFs, as well as Robin Stinchcombe and J. A. Castro for 
interesting discussions and suggestions. Thanks are also 
due to a referee for pointing out Ref. |23. This research 
was partially supported by the Brazilian agencies CNPq 
(Grant No. 30.0003/2003-0), FAPERJ (Grant No. E26- 
152.195/2002), FUJB-UFRJ and Instituto do Milenio de 
Nanociencias-CNPq. 



[1] A.-L. Barabasi and H. E. Stanley, Fractal Concepts 
in Surface Growth (Cambridge University Press, Cam- 
bridge, 1995). 

[2] M. Kardar, Phys. Rep. 301, 85 (1998). 

[3] S. T. Bramwell, P. C. W. Holdsworth, and J.-F. Pinton, 
Nature (London) 396, 552 (1998). 

[4] S. T. Bramwell, K. Christensen, J.-Y. Fortin, P. C. W. 
Holdsworth, H. J. Jensen, S. Lise, J. M. Lopez, 
M. Nicodemi, J.-F. Pinton, and M. Sellitto, Phys. Rev. 
Lett. 84, 3744 (2000). 

[5] T. Antal, M. Droz, G. Gyorgyi, and Z. Racz, Phys. Rev. 
Lett. 87, 240601 (2001). 

[6] T. Antal, M. Droz, G. Gyorgyi, and Z. Racz, Phys. Rev. 
E 65, 046140 (2002). 

[7] J. S. Urbach, R. C. Madison, and J. T. Markert, Phys. 
Rev. Lett. 75, 276 (1995). 

[8] M. Bahiana, B. Koiller, S. L. A. de Queiroz, J. C. 



Denardin, and R. L. Sommer, Phys. Rev. E 59, 3884 

(1999). 
[9] S. L. A. de Queiroz and M. Bahiana, Phys. Rev. E 64, 

066127 (2001). 
[10] S. L. A. de Queiroz, Phys. Rev. E 69, 026126 (2004). 
[11] H. Leschhorn, Physica A 195, 324 (1993). 
[12] H. A. Makse and L. A. N. Amaral, Europhys. Lett. 31, 

379 (1995). 
[13] H. A. Makse, S. Buldyrev, H. Leschhorn, and H. E. Stan- 
ley, Europhys. Lett. 41, 251 (1998). 
[14] A. Rosso, A. K. Hartmann, and W. Krauth, Phys. Rev. 

E 67, 021602 (2003). 
[15] G. Durin and S. Zapperi, e-print cond-mat 0404512. 
[16] E. Puppin, Phys. Rev. Lett. 84, 5415 (2000). 
[17] D.-H. Kim, S.-B. Choe, and S.-C. Shin, Phys. Rev. Lett. 

90, 087203 (2003). 
[18] K. L. Babcock and R. M. Westervelt, Phys. Rev. Lett. 



10 



64, 2168 (1990). 
[19] P. J. Cote and L. V. Meisel, Phys. Rev. Lett. 67, 1334 

(1991). 
[20] K. P. O'Brien and M. B. Weissman, Phys. Rev. E 50, 

3446 (1994). 
[21] S. Zapperi, P. Cizeau, G. Durin, and H. E. Stanley, Phys. 

Rev. B 58, 6353 (1998). 
[22] A. Rosso, W. Krauth, P. LeDoussal, J. Vannimenus, and 

K. J. Wiese, Phys. Rev. E 68, 036128 (2003). 
[23] G. Foltin, K. Oerding, Z. Racz, R. L. Workman, and 



R. K. P. Zia, Phys. Rev. E 50, R639 (1994). 
[24] T. Antal, M. Droz, and Z. Racz, J. Phys. A 37, 1465 

(2004). 
[25] K. Binder, Z. Phys. B 43, 119 (1981). 
[26] I. J. Zucker and M. M. Robertson, J. Phys. A 8, 874 

(1975). 
[27] S. Moulinet, A. Rosso, W. Krauth, and E. RoUey, Phys. 

Rev. E 69, 035103(R) (2004). 
[28] Z. Racz and M. Plischke, Phys. Rev. E 50, 3530 (1994). 



